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ABSTRACT 

We study the phase-space distribution of satelHte galaxies associated with late-type galaxies 
in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically rep- 
resentative regions from the Millennium simulation, which have higher resolution and incor- 
porate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions 
(r ~ 0.1r2oo) of the dark matter halo, both in shape and angular momentum, there can be 
substantial misalignments at larger radii (r ~ ?'20o)- Misalignments of > 45° are seen in 
^ 30% of our sample. We find that the satellite population aligns with the shape (and angular 
momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak 
owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite 
systems with ten bright galaxies within r-zoo exhibit a polar spatial alignment with respect to 
the galaxy — an orientation reminiscent of the classical satellites of the Milky Way. We find 
that a small fraction (~ 10%) of satellite systems show evidence for rotational support which 
we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the 
spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). 
This preference towards co-rotation is stronger in the inner regions of the halo where the most 
massive satellites accreted at relatively early times are located. 

We attribute the anisotropic spatial distribution and angular momentum bias of the satel- 
lites at z = to their directional accretion along the major axes of the dark matter halo. The 
satellite galaxies have been accreted relatively recently compared to the dark matter mass and 
have experienced less phase-mixing and relaxation — the memory of their accretion history 
can remain intact to z = 0. Understanding the phase-space distribution of the z = satellite 
population is key for studies that estimate the host halo mass from the line of sight velocities 
and projected positions of satellite galaxies. We quantify the effects of such systematics in 
estimates of the host halo mass from the satellite population. 

Key words: galaxies; general - galaxies; haloes - galaxies; kinematics and dynamics - dark 
matter - cosmology; theory 



1 INTRODUCTION 

Studies of local galaxies, such as the Milky Way and M31, can 
potentially provide us with the missing link between cosmological 
structure formation and the complex baryonic processes that help 
shape the galaxies we observe today. Any acceptable cosmological 
model must, in addition to satisfying the requirements of large scale 
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structure, account for the small scale detail exhibited by our own 
Milky Way galaxy and others. 

For some time it has been known that the eleven classical 
satellites of the Milky Way d efine a highly incl ined plane rela- 
tive to the disc o f the Galaxy (|L ynden-Bel| l 19761) . Early work by 
iHolmberd ( 1 19691) and IZaritskv et al. ( 19971) found similar align- 
ments in external galaxies whereby the satellites tend to avoid the 
equatorial regions of the parent light distribution. However, more 
recent work using the 2 degree field Galaxy Redshift Survey and 
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the Sloan Digital Sky Survey (SDSS) find that the opposite trend 
is true; satellites tend to avoid t he polar region s of the light distri - 
bution dSales & Lambasll2004 ISrainerdl l2005l : I Yang et alJlioO^ . 
lAgustsson & Brainerj ( I2OI0I) find that the spatial orientation of 
satellite galaxies depends on the type of host galaxy; satellites 
of red galaxies are found preferentially near the major axes of 
their hosts while satellites of blue galaxies show little or no spa- 
tial ansiotropy. The anisotropic nature of the spatial distribution 
of satellite galaxies has led some authors to postulate a discrep- 
ancy between observati onal constraints and models adopting the 
ACDM framework (e.g. iKroupa et al]|2005h . However, numerical 
studies have shown that p referential alignments are naturally pro- 
duced in the simulations l iKang et al.||20o3; iLibeskind et ai]|2005l ; 
IZentner et al .l2005h . The satellites tend to align with the major axis 
of the dark matter halo, but an extrapolation to the relation with 
the light distribution is not straightforward in models which do not 
follow the evolution of baryonic matter. 

The radial velocities and available proper motions of the clas- 
sical Milky Way satellites hint ait the presence of cohe rent mo- 
tion llLvnden-Bell & Lvnden-Bell|[T99l ; lMetz et alj|2008l) . In fact, 
iMetz et alj j2008t) find that at least 3 of the classical satellites have 
orbital poles aligned (within 30°) with the normal of their spatially 
defined plane. They suggested that the classical satelli tes may oc- 
cupy a rotationally-supported disc. iLi & Helmil ( |2008|) found that 
coherence in the motions of the satellites may be due to group in- 
fall, whereby satellites fall into the parent halo together and pre- 
serve their common motion to the present day. A scenario whereby 
satellite accretion is along surrounding filamentary structures sug- 
gests a link between the angular moment um orientation of the 
satellite galaxies and the host halo itself ( lLibeskinderai]|2005l ; 
iLovell et al.ll20T l|). A correlation between the orbital motion of 
satellit e galaxies and the sp in of their parent light distribution was 
seen bv lAzzaro et al.l ( l2006h in a carefully selected sample of SDSS 
galaxies. However, th is bias towards co-rotation was not seen by 
iHwang & Pa^ ( I2OI0I) in a larger sample of SDSS host galaxies. 
The authors find equal numbers of satellites in prograde and ret- 
rograd e or bits. Numeri c al sim ulation studies bv lWarnick & Knebg 
( I2OO6I) and lShaw et all ( l2006h find a bias towards satellites on co- 
rotating orbits relative to the net spin of their host halo. How- 
ever, both of these studies base their conclusions on results derived 
from dark matter only simulations, and they focus on cluster-sized 
haloes. Observational estimates are frustrated by uncertainties re- 
garding the spin direction of the parent galaxies and contamina- 
tion by interlopers. On the other hand, theoretical work has solely 
focused on dark matter only simulations — the orientation of the 
satellites' orbit with respect to the stellar distribution is yet to be 
tested. 

Dark matter only simulations have been hugely influential in 
developing our knowledge of the large scale structure of the Uni- 
verse. However, some of the potential shortcomings of the standard 
model posed by observations of our own Milky Way galaxy are dif- 
ficult to reconcile within a simulation that does not include lumi- 
nous matter. Previous work has made use of semi-analytic models 
in order to include the necessary b aryonic processes into these cos- 
mological simulations (e.g. White & Frenk|[T99ll ; ICole et al.lll994l ; 
ISomerville & Primack|[T999, ; .Baughll20oj) . Whilst these methods 
have provided valuable insights into the effects of baryonic physics 
on galaxy formation, their limited spatial information makes them 
unsuitable to investigate the processes of interest in this paper. For 
example, a key assumption in these semi-analytic methods is that 
the angular momentum of the disc is aligned with the spin of the 
dark matter halo. 



In recent years, increasingly realistic implementations of 
the hydrodynamic evolution of the baryons have become possi- 
ble within the framework of cosmological simulations which are 
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lar, these hydrodynamic cosmological simulations are able to fol- 
low the changes in shapes and angular momenta of the dark mat- 
ter and baryons self-consistently. In this paper, we make use of 
the GI MIC suite of simulations described in detail bv ICrain et al.l 
( l2009h . This is a re-simulation of 5 cosmologically representative 
regions (^ 20fe~^ M pc in radius) from the Millennium simulation 
jSpringel et allbOOSh . In these regions, GIMIC incorporates bary- 
onic physics and achieves higher resolution than the Millennium 
simulation as a whole. 

We use these simulations to study the orbital properties of 
satellite galaxies in late-type galaxies. In contrast to previous work, 
we probe the dynamics of the satellites relative to their host's stel- 
lar component as well as the unseen dark matter component. In §2, 
we describe the GIMIC suite of simulations in more detail and out- 
line the selection criteria for our sample of parent haloes. In §3, we 
discuss the relation between the galaxy and the dark matter halo in 
both spatial and velocity space. §4 focuses on the satellite galaxies 
associated with our sample of haloes. We investigate their spatial 
and angular momentum distribution relative to both the galaxy and 
dark matter halo. In §5, we briefly outline an application of our re- 
sults to test estimators of the parent halo mass. Finally, in §6 we 
draw our main conclusions. 



2 THE NUMERICAL SIMULATIONS 

In this section, we briefly describe the simulations we have used 
and outline our methods for selecting parent galaxy haloes and their 
associated satellite galaxies. 



2.1 GIMIC 

The Galaxies-Intergalactic Medium Interaction Calculation 
(GIMIC) suite of simulations is des cribed in detail in ICrain et al.l 
l l2009h (see also ICrain et al.l I2OI0I) . It consists of a set of hy- 
drodynamical re-simulations of five nearly spherical regions 
(~ 20fe ~^ Mpc in radius) extracted from the Millennium Sim- 
ulation jSpringel etal.ll2005l . The regions were selected to have 
overdensities at z — 1.5 that represent (+2, +1, 0, — 1, — 2)a, 
where a is the root-mean-square deviation from the mean on 
this spatial scale. The 5 spheres therefore encompass a wide 
range of large-scale environments. In the present study, we select 
systems with total 'main halo' (i.e., the dominant subhalo in a 
friends-of-friends group) masses s imilar to tha t of the Milky Way, 
irrespective of the environment. ICrain et aL I l l2009l) found that 
the properties of systems of fixed main halo mass do not depend 
significantly on the large scale environment (see, e.g.. Fig. 8 of that 
paper). 

We present only a brief summary of the GIMIC simulations 
here, and refer to Grain et al. (2009, 2010) for more detailed de- 
scriptions. The cosmological parameters are the same as those in 
the Millennium Simulation and correspond to a ACDM model with 

= 0.25, = 0.75, f7b = 0.045, as = 0.9 (where erg is 
the rms amplitude of linear mass fluctuations on 8h~^ Mpc scale 
at z = 0), Ho = lOO/i km s"^ Mpc"\ h = 0.73, Us = 1 (where 
Us is the spectral index of the primordial power spectrum). 
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Table 1. The average properties of our sample of (431) galaxy haloes. We 
give the median halo radius (r200 ), the median total mass inside this radius 
(A^20o). the median parent halo stellar mass (Mi,), the median disc-to-total 
stellar mass ratio (D/T) and the median number of satellite galaxies per 
parent halo (A^sat)- 



The simulations were evol ved to z = u sing the TreePM- 
SPH code GADGET, described in lSpringeJllOOSl . Subsequently, the 
GADGET code has been substantially modified to incorporate bary- 
onic physics which includes: 

• a prescription for star formation outlined in 
ISchave & Dalla Vecchij jlOOSh tha t is designed t o reproduce 
the observed Kennicutt-Schmidt law l lKennicutti[l998l) : 

• radiative gas cooling in the pres ence of a UV/X-Ray back- 
ground (see iHaardt & Madaul 1200 ll) which includes the con- 
tribution of metal-lin e cooling (computed element-by element) 
dWiersma et alj2009ah : 

• the timed release of 11 individual metals by both massive 
(Type II SNe and stellar winds) and intermediate mass stars (Type 
lA SNE and asymptotic giant branch stars ) ( Wiersma et al. 20 09bO; 

• a kinetic supernova feedback model jPalla Vea:hia & Schavd 
l2008h which can quench star-formation in low mass haloes and pol- 
lute the IGM with metals. 

In the present study we analyse the 'intermediate' resolution 
GIMIC simulations which have 8 times better mass resolution than 
the Millennium simulation. These runs have a dark matter parti- 
cle mass Mdm ~ 5.30 x lQ^h~^ M© and an initial gas particle 
mass of Mg = 1.16 x 10'^h~^ M©. This implies that it is possi- 
ble to resolve systems with masses similar to that of the classical 
dwarf galaxies. The remainder of the 500(/i~^ Mpc)"^ Millennium 
volume is modelled with dark matter particles at much lower reso- 
lution to ensure the presence of surrounding large scale structure is 
accurately accounted for. Dark matter only runs of the GIMIC sim- 
ulations are also available (at the same resolution) which can be di- 
rectly compared to the hydrodynamic versions of the simulations. 
The — 2cr volume of the 'high' resolution GIMIC simulations has 
also been carried out to redshift z = Q. This higher resolution run 
has 8 times better mass resolution than the intermediate runs (and 
64 times better than the Millennium simulation). We check that the 
conclusions of this paper are not subject to resolution effects by 
ensuring our main results are unchanged in the higher resolution 
simulations (see Appendix A). 

Previous work utilising the GIMIC simulations has found en- 
couraging agreement with observational studies. The adopted star 
formation and feedback prescription results in a good match to 
the star formation rate histo ry of the universe (Grain et al. 2009; 
see also ISchaveetal.ll2010l) . and also reproduces a nu mber of X- 
ray/o p tical scaling relat ions for normal disc galaxies jCrain et al] 
|201(]|) .I Font et al. find that the stellar haloes of their simu- 

lated Milky way-mass galaxies have luminosities and radial density 
profiles in good agreement with observations. Furthermore, Mc- 
Carthy et al. (In prep) show that the simulated L* disc galaxies in 
GIMIC have realistic kinematics and sizes. 



2.2 Identification of Galaxies and Satellites 

Bound haloes are identified using the SUBFIND algorithm of 
Dolag et al.l ||2009|), w hich extends the standard implementation of 
Springel et al. I l l200lh 3y also including baryonic particles when 
identifying self-bound substructures. The main galaxy or parent 
halo is the most massive subhalo belonging to a friends-of-friends 
system. The other self-bound substructures of the system are then 
classified as satellite galaxies. 

We define r2oo as the radius at which the mean enclosed 
density falls to 200 times the critical density (200pcrit)- We se- 
lect parent haloes with total mass within this radius in the range 
5 X 10" < A/200 /M© < 5 X 10^^. Only haloes with at least 
one associated satellite galaxy (or subhalo) other than the main 
halo are included. GIMIC is a re-simulation and is thus subject to 
edge effects at the boundaries of the selected spherical regions. We 
discard any haloes located at the boundary edges that are partially 
comprised of low resolution dark matter particles. 

We select our sample of parent haloes according to the 're- 
laxation' criteria defined below, following the same reasoning as 
iNetoetallfcoOTl) : 

(1) Virial ratio: We compute the ratio 2K/\U\, where K and U are 
the total kinetic energy and total potential energy within r2oo- We 
adopt 2J\/|(7| < 1.35 for the relaxed sample. 

(2) Centre of mass displacement: The offset between the centre of 
mass of the halo and position of the most bound particle (potential 
centre) can be de scribed by the norm alised offset parameter, e = 
\tc — rcm\/r200 l lThomas et al.ll20()l]) . Relaxed haloes have e < 
0.07. 

(3) Mass in substructure: For a halo to be considered relaxed, 
we require the fraction of mass in substructure within r2oo to be 
/sub < 0.1 and the most massive satellite within 2r2oo to satisfy 
Msat/Maoo < 0.1. 

Under t hese criteria appro ximately 20% of the haloes are un- 
relaxed (c.f. iNeto et ai]|2007h . There are a total of 624 relaxed 
haloes within the given mass range. Most unrelaxed haloes are 
recognised from the second and third constraints, whilst only a few 
haloes are out of virial equilibrium. Note that since only a rela- 
tively small fraction of haloes are unrelaxed, our conclusions are 
generally valid. 

The galaxies are assigned a morphological classification (i.e. 
disc and spheroid-dominated types) based on their dynamics. 
A simple two-component model is assumed: (i) a dispersion- 
supported spheroid, and (ii) a rotationally supported dis c. Details of 
this de compositio n into morphologic al types is given in lCrain et al.l 
( I2OIOI) . Fig. 1 in lCrain et al.l j2010l) shows the distribution of the 
disc-to-total stellar mass ratios (D/T). In this work, we restrict 
our sample to late-type galaxies and exclude obviously 'elliptical' 
galaxies. Following the reasoning of ICrain et alj fcOlOl) . we adopt 
a threshold of D/T > 0.3. We have checked that our main results 
are unchanged if other cuts of D/T are imposed (e.g. D/T > 0.2 
or D/T > 0.4). We note that haloes are selected based only on 
mass and D/T. This selection was not chosen to reproduce Milky 
Way (or M3 1 ) galaxy replicas and comparison with these two Local 
Group galaxies is made in the broadest sense. 

Our final sample of haloes consists of 43 1 parent haloes and 
4864 associated satellite galaxies. We summarise the properties of 
our sample in Table [T] The approximate (total) mass range of our 
sample of satellite galaxies is W^Mq < M^at < lO^Af©. Pre- 
vious work using low resolution dark matter only simulations have 
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Figure 1. The radial dependence of the axial ratios s = c/a (top panels), p = c/b (second row), q = b/a (third row) and the triaxiality parameter 
T = (1 — g^)/(l — s^) (bottom panels). The shaded regions show the range of values covered by 68% of the distribution and the points give the median 
values. The dotted lines indicate the regions of oblate (T < 1/3), triaxial (1/3 < T < 2/3) and prolate haloes (T > 2/3). The right hand panels are for the 
dark matter only counterparts of our sample. The halo shapes are rounder (and more oblate) in the hydro GIMIC simulations. 



much larger samples of parent haloes (e.g. ISales et al.ll2007L who 
use thousands of parent haloes from the Millennium simulation). 
However, our sample size compares favourably to more recent hy- 
drodynamical simulations (e.g. Libeskind et al. 2007, who had only 
9 parent haloes). 

In addition to our hydrodynamical suite of GIMIC simulations, 
we have also the dark matter only runs. We can match our sample of 
haloes with their dark matter only counterpart. For each halo, we 
consider dark matter only haloes with similar masses inside r2oo 
(within a factor of 2). We then compute the distances between the 
positions of the dark matter only haloes and the position of the halo 
in question. The dark matter only counterpart is then the closest (in 
position) to the baryonic simulations version of the halo which also 
has a similar mass. 



(or reduced) second moment tensor 

]V(s;r) 



where 



+ 



''2.n 



+ 



''3,n 



(1) 



(2) 



The advantage of this scheme is that every particle is given equal 
weight, independent of radius. The orientation is defined by the 
eigenvectors of the second moment tensor. The principal compo- 
nents of the tensor give the axial ratios 



(3) 



3 HALO PROPERTIES 

Here, we examine the shapes, the angular momenta and the mis- 
alignments of the dark matter haloes of our simulated galaxies. 



3.1 Shapes 

Haloes are modelled as ellipsoids characterised by three axes, 
a,b,c where a ^ b ^ c. The axial ratios s = c/a, q — b/a 
and p = c/b describe the three classes of ellipsoids: prolate 
(a > b ~ c), oblate (a ~ b > c) and triaxial (a > b > c). 
The shapes of dark matter haloes are computed using the weighted 



As the value of the elliptical radius, r„, is not known in advance (it 
depends o n s and q), the axial ratios are computed using an iterative 
algorithm jPubinski & Carlberg|[l99ll) . In Figure [T] we show the 
radial dependence of the three axis ratios. The bottom panel gives 
the radial behaviour of the triaxiality parameter ( lFranxetal.lll99lh 
defined by 



[b/af 



{c/aY 



(4) 



Oblate, triaxial and prolate haloes have triaxiality parameters of 
r < 1/3, 1/3 < T < 2/3 and r > 2/3 respectively. The right 
hand panels of Fig. [T] give the axis ratios for the dark matter only 
counterparts of our sample. The shaded regions show the values 
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Figure 2. Median misalignment angles as a function of radius. EiTor bai's and shaded regions sliow the values covered by 68% of the distribution (i.e. the Icr 
dispersion). Top left: Misalignment angles between the angular momentum vector of the galaxy (jgai) and the angular momentum vector of the dark matter 
halo (jdm)- Top light: Misalignment angle between the short axis of the galaxy (Cg^i) and the short axis of the dark matter distribution (cj„j). Bottom left: 
Misalignment angles between the major axis of the galaxy (Sgai) and the major axis of the dark matter distribution (Sdm)- Bottom right: MisaUgnment angles 
between the short axis and the angular momentum vector of the dark matter distribution (green shaded region). For comparison the misalignment between the 
short axis of the velocity anisotropy tensor (Cvd) and the short axis of the dark matter halo shape (Cghapc) is shown by the red shaded region. There are strong 
alignments between the galaxy and the dark matter halo at small radii (r ~ 0.1r20o) but there can be significant misalignments at larger radii (r ~ ^200)- 



covered by 68% of the distribution and the points denote the median 
values. 

The dark matter haloes in our sample are close to spherical and 
are slightly more oblate in the inner regions. This is in stark contrast 
to the haloes in the dark matter only simulations which become 
more prolate towards the centre. The inclusion of baryonic physics 
affects the shape of the halo substantially in the inner regions, but 
also has a significant effect throughout the halo. 

Our results broadly agree with the findings of previous studies 
but direct comparisons are difficult as many authors have concen- 
trated on cluster-sized haloes or have used simulations where the 
assembly of the central galaxy is not modelled self-consistently. 
iDubinskil l ll994) . for example, simulates dissipative infall by grow- 
ing a central mass concentration inside a triaxial dark matter halo 
and find that a steeper potential leads to rounder and more oblate 
dark matter halo shapes (A(c/a) ~ O.D. lKazantzidis et alj ^20041) 
study cluster-sized haloes using gas dynamical simulations and find 
they are significantly rounder in the inner regions (A(c/a) ~ 0.2) 
but the changes are ra dially depen d ent an d are almost negligi- 
ble at the virial radius. lAbadi et"al] 1 I2OIOI) employ galaxy mod- 
els which include radiative gas cooling but neglect the contribution 
from stellar feedback. They compute equipotential axial ratios and 
find roughly constant flattening of (c/a) ~ 0.85, which is signif- 
icantly rounder than their dark matter only runs. By contrast, we 
have characterised the halo shape by the density of the dark matter, 
w hich is always flatt er than the equipotential surfaces. Compared 
to lAbadieta"n i l2010l) . we do not find such a significant change in 



sphericity relative to the dark matter only simulations. This is pre- 
sumably because the GIMIC simulations include stell ar feedback 
and d o not suffer from strong overcooling (see e.g. iDuffv et al] 
|201(]|) . As indicated in Table [T] the mean stellar mass fraction of 
our simulated disc galaxies is ~ 0.057. This corresponds to a 
baryon conversion efficiency (M,/A/haio x Qm/^b) of « 30%, 
which is only a f actor of ~ 1.5 larger than that inferred recently by 
IGuo et al.l ( I2OI0I) by matching the observed stellar mass function 
from the most recent SDSS data release to the halo mass function 
derived from the Millennium and Millenniwn-II simulations. 

Unfortunately, the evidence on the shape of the Milky Way's 
dark halo is far from clear-cut. Many authors have argued that 
the coherence of the Sagittarius tidal stream may constrain the 
halo shape. However, this line of enquiry has co n cluded that the 
halo may be almost spherica l 1 Ibata et aP '2OOI'; Fellhau er et al.l 
l200db. oblate jjohnstoii et al.ll2005h . prolate (Helmi 2004i) or tri- 
axial jLaw & Maiewskill2010h . This variety of results strongly sug- 
gests that halo shape is not the primary facto r determining the c om- 
plex morphology of the Sagittarius stream. Ismith et afl ( |2009|) ar- 
gued that the spherical alignment of the velocity ellipsoid of SDSS 
halo subdwarf stars implied that the gravitational potential and 
hence the dark halo was nearly spherical. There is however con- 
tra dictory evidence from stud ies of the flaring of the HI gas layer 
bv lOlling & Merrifielddioool) . who found a highly flattened oblate 
(g ^ 0.3) dark halo for the Milky Way. 

The dark matter halo shapes of external galaxi es are often es- 
timated using galaxy-galaxy weak lensing studies. iHoekstra et al.l 
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Figure 3. Top: The distribution of misalignment angles between the angular momentum of the galaxy and the angular momentum of the dark matter halo for 
r ^ 0.1r200 (left panel) and r r200 (right panel) respectively. Bottom: The distribution of misalignment angles between the short axis of the galaxy and 
the short axis of the dark matter halo for r 0.1r2oo (left panel) and r r2oo (right panel) respectively. The dotted lines show a uniform distribution. 



( I2OO4I) and lParker et"aLl j2007h both found an average axis ratio of 
~ 0.7 from their studies using the Red-Sequenc e cluster survey and 
the CF HT legacy survey, respectively. However. lMandelbaum et all 
( I2OO6I) found no evidence for halo ellipticity in their study using 
SDSS data. 



3.2 Angular Momentum and Shape 

The (cumulative) specific angular momentum vector for the dark 
matter particles, jdm, and the stellar component, jgai, are defined 
as 

^ iV(^r) 

»•) ^ MK r) ^ "i„x„xv„ (5) 

n — l 

where Xn and Vn are the position and velocity vectors of parti- 
cle n relative to the halo centre and the centre-of-mass velocity. We 
compute jgai for r ^ 0. lr2oo (using only the star particles) to char- 
acterise the inner galaxy. The stress or velocity dispersion tensor of 
the dark matter distribution is 

Ilij = ^ Vi^nVj,n- (6) 
n=l 

We can diagonalise this tensor to find the eigenvectors and 
eigenvalues which define the principal velocity anisotropy axes 

(Ct vel , ^vel 7 Cvel ) . 

In Fig. |2l we show the median misalignment angles between 
the galaxy and the dark matter halo as a function of radius. The 



top left panel shows the misaligmnent in angular momentum, the 
top right the misalignment in shape. Note that the angular momen- 
tum or shape of the galaxy is always calculated for r ^ 0.1r2oo, 
whereas the computation of the dark matter halo properties varies 
with radius. There is strong alignment of angular momentum vec- 
tors in the inner regions of the halo, but the median misalignment 
gradually i ncreases to ~ 40 ° at r ~ r2oo. This is in good agree- 
ment with iBett et al] l l201(]h who find a median misalignment of 
~ 30° for a smaller sample of ~ 50 haloes. There is also strong 
alignment between the short axis of the galaxy and the short axis 
of the dark matter halo where the median misalignment grows to 
~ 20° at r ~ r2oo. In the bottom left hand panel, we see that the 
alignment between the major axes of the galaxy and dark matter 
halo is poorer and there is significant scatter. This is not surprising 
as for disc-like configurations the major axis is poorly defined (i.e. 
q = h/a ~ 1). The bottom right hand panel shows the alignment 
between the short axis and the angular momentum vector of the 
dark matter halo as a function of radius (shaded green). We also 
show as the shaded red region the alignment between the short axis 
of the dark matter halo shape and the short axis of the velocity dis- 
persion tensor. There is very strong alignment between the shape 
and velocity dispersion. This reflects the fact that the dark matter 
halo shape is supporte d by internal velocities rather than net rota- 
tion jprenket all 19881) . 

To avoid any ambiguity in the definition of the sh ort axes or 
angula r momentum vectors of the dark matter haloes, iBett et all 
l l201(]h imposed constraints on their shapes and net angular mo- 
mentum (s < 0.81, logioj(^ J"2oo)/VGM2oo'"2oo 5= -1.44). 
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With these restrictions, we see stronger alignments in general be- 
tween the galaxy and dark matter halo. This is not surprising, as 
we have found that the dark matter haloes are close to spherical 
and are mainly dispersion, rather than rotationally, supported. How- 
ever, we choose to include all of our haloes to avoid selection bi- 
ases (only ~ 40% of our sample satisfy the constraints), but check 
that our main results are not significantly altered by discarding dark 
matter haloes that have little net rotation and are close to spherical. 

Further insight into the misalignment between the galaxy and 
dark matter halo can be gained by examining the distribution of 
these misalignment angles at different radii. In Fig. [3] we show the 
distributions of misalignments between the galaxy and dark matter 
halo angular momentum vectors (top panels) and short axes (bot- 
tom panels) both for r ^ 0.1r2oo (left panels) and for r ^ r2oo 
(right panels). As we already saw in Fig.|2] there is much stronger 
alignment in the inner regions of the halo. However, when the shape 
of the dark matter halo is computed for r ^ r2oo, the short axis 
of the inner galaxy is misaligned from the short axis of the dark 
matter halo hy 6 > 45° in approximately 30% of the systems. 
The orientation of the angular momentum vector of the galaxy is 
aligned almost perpendicularly (45° < 9 < 135°) to the net spin 
of the dark matter halo in approximately 40% of the systems and 
2% are even anti-al igned. Our results are in good agreement with 
iBailin et al. I l l2005al) who analyze the shape and internal alignment 
of the dark matter halos of seven high-resolution cosmological disk 
galaxy formation simulations. The authors find that while the mi- 
nor axis of the inner dark matter halo (r ^ 0.1r2oo) is well aligned 
with the disc axis, the outer regions of the halo can be substantially 
misaligned. 

The misalignments between the galaxy and the dark matter 
halo have potential implications for galaxy-galaxy weak lensing 
studies. Such studies attempt to deduce the weak cosmological sig- 
nal from the observed galaxy ellipticity correlation. Typically the 
influence of any intrinsic signal on the observed ellipticity corre- 
lation is neglected. However, correlations between halo shape and 
the density field are expected to arise through tidal torques operat- 
ing during galaxy formation (e.g. lHeavens & Peacock! 19881) . More 
recent work has attempted to quantify the potential sources of lens- 
ing contamination caused by these intrinsic alignm ents by examin- 
ing d ark matter only cosmological simulations (e.g. lHevmans et all 
l2006h . Any misalignments between the galaxy and the dark matter 
halo will weaken the expected intrinsic signal deduced from dark 
matter only simulations. 

What causes these misalignments between the galaxy and 
the outer dark matter halo? According to tidal torque theory (e.g. 
[white 1984; Fall & Efstathiou 1980), the angular momentum of 
the galaxy and the dark matter halo are initially very well aligned. 
However, the outer halo continues to accrete material, which can 
alter its shape and/or its net angular momentum. Hence, whilst we 
see strong alignment in the inner regions of the halo there can be 
significant misalignments in the outer parts. This has important im- 
plications for the satellite populations of these haloes. The majority 
of the satellite galaxies are located in the outer regions of the halo. 
It is to this topic that we now turn. 



4 SATELLITE GALAXIES 

Here, we study the spatial and velocity distributions of the satellite 
systems of our simulated galaxies, focusing on their alignment (in 
positional and velocity space) with the parent dark matter halo and 
galaxy. 



satellites: a = 0.62, r^ = 0.40 
dark matter; a = 0.16, r =1e-6 
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Figure 4. Top panel: The cumulative number distribution for satellites 
(full red line) and dark matter particles (full blue line) within r200- Fits 
to Einasto profiles are given by the red and blue dashed lines. The dot- 
ted magenta and cyan lines show the cumulative number distribution for 
satellites with stellai' mass fractions, M*^^ / Af sat , greater than 2% and less 
than 2% respectively. Middle panel: The red and blue lines give the radial 
profiles of the satellites and dai'k matter particles respectively. The dashed 
Hues show the profiles for the dark matter only counterparts of our sample 
of haloes. Bottom panel: The cumulative number distribution for the ten 
brightest satellites of individual systems. The blue hashed region encom- 
passes the scatter of the profiles for the ~ 80 systems with ten or more 
satellites within r200 ■ The red and green lines give the profiles for the clas- 
sical Milky Way satellites and the ten brightest M31 satellites (within r2oo) 
respectively. 
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In Figure |4] we show the cumulative number distribution of 
satellite galaxies (red) and dark matter (blue) within r2oo- We have 
stacked all of the satellite galaxies by normalising their radial dis- 
tances from the parent halo centre by r2oo • The density profiles for 
both the dark matter particles and satellite galaxies can be described 
jNavarro et"ar]|2004) by 

Inp(r) oc -(2/Q)[(r/r,)" - 1]. (7) 

This density profile was first introduced by Einasto (see e.g. 
lEinasto&Haiidlll989h and is mathematically equivalent to the Ser- 
sic profile that is often used to describe the projected density pro- 
file of galaxies. Low values of the power-law slope (a) describe 
a 'cuspy' density profile. In agreem ent with previous work (e.g. 
iNavarro et aLll2004 lOao et allboOSh , we find a dark matter den- 
sity slope of Qdm ~ 0.2. In contrast, the satellite distribution 
favours larger values o f a and has a 'cored' profile. As found 
many time s before (e.g. lZentner et al ] |2005l : iLibeskind et allbOOSl : 
iLudlow et al. 2009), there is an obvious spatial bias between the 
dark matter and the satellite galaxy population. 

Fig.|4]also shows the cumulative number of satellites with stel- 
lar mass fractions greater than 2% and less than 2% by the ma- 
genta and cyan dotted lines respectiveljQ . Satellites with higher 
stellar mass fractions are more centrally located. These are gener- 
ally the most massive satellites which have spiralled into the central 
regions of the halo during the course of several pericentric pas- 
sages. The stellar component of satellite galaxies is more centrally 
located than the dark matter distribution. Whilst the dark matter 
'envelope' is more effectively tida lly stripped, the stellar compo- 
nent remains relatively undisturbed. ILibeskind et"aLl ( 1201 oh showed 
that the increased density in the central regions of a subhalo due to 
the collapse of baryons leads to a reduced m ass loss relative to a n 
equivalent dark matter only satellite (see also lMaccio et aklboOfil) . 
Dynamical friction is thus more effective and satellites sink closer 
into the parent halo. This leads to a more centrally concentrated dis- 
tribution of satellite galaxies when baryons are included. By com- 
parison with our dark matter only simulations, we also see a more 
substantial inner radial bias when baryons are included (see bottom 
panel of Fig. |4j. Note that this effect is not as pronounced for the 
lower mass satellites, which are less affected by dynamical friction. 

The bottom panel of Fig. |4] shows the cumulative number dis- 
tribution of the ten brightest (i.e. highest stellar-mass) satellites 
within r2oo for individual systems. Approximately 80 haloes have 
at least ten satellites within 7-200 ■ We limit ourselves to the ten 
brightest satellites to enable comparison with the data on the Milky 
Way and M3 1 . The positional data for the Milky Way and M3 1 
satellites is taken from tables 3 and 5 of lOeason et alj ( I2OI l[) (see 
references therein). Although the top panel of Fig|4| gives no indi- 
cation of the scatter in the simulations, the bottom panel explicitly 
shows the system-to-system variation (encompassed by the blue 
hatched region). The solid red and green lines gives the profiles for 
the classical Milky Way satellites and the ten brightest M31 satel- 
lites (within r-2oo)- We use the r2oo values f or the Milky Way and 
M31 recently estimated bv lGuo et al. I ( I2OIOI) of approximately 250 
kpc and 290 kpc respectivel50. The profiles for both M3 1 and the 
Milky Way lie within the scatter of the simulations. In particular, 
there is very good agreement with M3 1 for which (owing to our ex- 



^ We verified that our results are not substantially affected if we choose 
fractional stellar mass limits of 1.5 or 2.5% 

^ These were comp uted using the estimated halo masses given in Table 1 . 
of lGuo et al.l j2010l) and then their eqn (1) was used to estimate r200 



temal view) there are less biases in the observed sample. Note that 
iFont et al] 1 I2OIII) show in their Fig. 1 that the luminosity function 
of the satellite galaxies in GIMIC are also in good agreement with 
the satellite systems of Local Group galaxies. 



4.1 Anisotropic Distributions 

Here, we first restrict attention to host galaxies which have ten 
or more associated satellites within 7-200 (cf. Libeskind et al. 
2007, 2009). For comparison with the classical satellites of the 
Milky Way (and M31), we consider the ten brightest (i.e. high- 
est stellar mass) satellites in each system. We compute the shape 
(flsat , &sat , Csat ) of the Satellite galaxy distribution within r2oo of 
each host galaxy by diagonalising the second moment tensor de- 
fined as 

JV{^r) 

^ij ~ ^ ^ ^i,n^j,n- (8) 
n=l 

This is used in preference to the reduced inertia of eqn ([TJ which 
requires an iterative algorithm to discard outliers until convergence 
is achieved. For systems of satellites with a small number of data 
points, this is undesirable. 

In the left hand panel of Fig.|5] we show that the satellite dis- 
tribution is generally more flattened than the underlying dark mat- 
ter distribution. The red dot indicates the median axial ratio values 
for the satellites and dark matter ((ssat) ~ 0.5, (sdm) ~ 0.8). We 
show with the blue arrow the range of values for the flattening of the 
Milky Way dark halo given in the literature. We suggest the differ- 
ence in flattening may be understood by considering when the dark 
matter was accreted relative to the satellite galaxies. Present-day 
satellites are the surviving population, and have have been in orbit 
for much less time than the dark matter (see Section |4!4l l. As far 
as the satellites are concerned, the potential of the halo has largely 
been static since accretion. The dark matter, by contrast, has under- 
gone relaxation and phase mixing. Even though both the satellites 
and dark matter are accreted anisotropically, this will be reflected 
to a greater degree in the spatial distribution of the satellite galax- 
ies rather than the dark matter halo. The right hand panel of Fig.|5] 
shows the distribution of axial ratios for the satellite systems. The 
blue and green stars show the flattening for both the classical Milky 
Way satellites and the ten brightest M3 1 satellites computed in the 
same way as the simulated satellites. The classical Milky Way satel- 
lites have a highly flattened distribution which is consistent (within 
2cr) with the median value derived from the simulations whilst the 
M31 satellites have a less flattened configuration and lie closer to 
the median s — c/a value. The points are coloured according to the 
alignment of the satellite distribution with the galaxy (see below 
and Fig.le}. Polar (Cgai.Csat < 0.25), planar (cgai.Csat > 0.75) and 
'in-between' (0.25 < Cgai.Csat < 0.75) alignments are given by 
the magenta, cyan and black points respectively. There is no obvi- 
ous correlation between the orientation of the satellite distribution 
and their flattening. Neither does there seem to be a bias towards a 
particular orientation, as we discuss below. 

We show the distribution of alignments between the short axes 
of the satellite systems and the short axes of the dark matter halo 
(for r ^ r2oo) and the galaxy of their parent haloes in Fig. |6] The 
orientation of each satellite system with respect to its host dark mat- 
ter distribution and galaxy are shown in the right and left panels re- 
spectively. The satellite distribution preferentially aligns in a plane 
perpendicular to the short axis of the dark matter distribution. How- 
ever, the satellites show no preferential alignment relative to the 
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Figure 5. Left panel: The shape of the satellite distribution and underlying dark matter haloes as defined by the sphericity, s = c/a. Only haloes with ten or 
more satellites within r20() are considered and the shape is computed for the ten highest stellar mass satellites. The red dot indicates the median values whilst 
the blue line indicates the approximate values for the Milky Way (the blue an'ows denote the uncertainty in the sphericity of the dai'k matter halo of the Milky 
Way). The points are colour coded according to the alignment between the system of satellites and the galaxy. Cyan, magenta and black points are for satellite 
systems where Cgai-Csat > 0.75, Cgai-Csat < 0.25 and 0.25 < Cg^i-Csat < 0-75, respectively. The satellite galaxies generally have a more flattened spatial 
distribution than the dark matter. Right panel: The axial ratios of the satellite distribution. The blue and green stars indicate the axial ratios for the classical 
Milky Way satellites and the ten brightest M31 satellites (within r20o)- 




Figure 6. The orientation of the short axes of the satellite distribution relative to the short axes of the parent galaxy and dark matter halo (defined within r2oo)- 
The black lines are for all systems of satellites with ten or more satellites within r200 . The red dashed lines give the distributions when only the ten brightest 
satellites within r2oo are used to compute the shapes of the satellite distribution. Downward pointing aiTows denote the median of the distributions and the 
dotted lines indicate a uniform distribution. The error bars denote Poisson uncertainties. Whilst there is strong alignment between the spatial distribution of 
the satellite galaxies and the dark matter halo, the distribution with respect to the galaxy is consistent with isotropy. 



galaxy. There are a significant number of systems where the satel- 
lite distribution is aligned in a plane perpendicular to the disc (20% 
of the systems have cos^ < 0.2 or 6 > 80°). Thus the alignment 
of the Milky Way satellites perpendicular to the disc plane is not 
inconsistent with the results we present here. The red dashed lines 
give the distributions when only the ten brightest satellites within 
r2oo are used to compute the shapes of the satellite distribution. 
Although there is substantial uncertainty caused by small number 
statistics, we can see that by restricting ourselves to the same num- 
ber of satellites as the observational sample (and naively ignoring 
selection biases), the apparent distribution of Milky Way satellites 
is consistent with the simulations. 

In Table[2l we give the number of satellite systems which have 
similar (perpendicular) orientations as the Milky Way Galaxy. Our 
results are robust to more restrictive D/T cuts. In fact, the frac- 



tion of systems with almost perpendicular satellite alignments is 
slightly higher for parent galaxies with D/T > 0.7. Our total mass 
range (5 x 10^^ < M-zoo/Mq < 5 x 10^^) is broadly coincident 
with the masses of the Milky Way and M31. We can further refine 
our sample selection by using haloes with stellar masses similar 
to these local galaxies. In the bottom row of Table |2l we give the 
number of satellite systems with polar alignments which also have 
parent stellar masses within a factor of 2 of the observed estimates 
for theMilky Way andM31 (2.5 x < M*/Mp, < 2 x lO" ; 
e.g. lWidrow et alJl2003l : rGeehan et al.ll2006l : lHammer et alj|2007h . 
Independently of our D/T cut, we find that approximately 20% of 
disc galaxies with similar masses as the Milky Way and M3 1 have 
satellite systems with polar alignments relative to the galaxy. 

In Fig.|7]we show the probability distribution of the orienta- 
tion of the satellite galaxies relative to their host. This differs from 
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Table 2. The number of parent haloes in the mass range 5 X 10^^ < 
M200 < 5 X 10^2 for different D/T cuts. We give the total number of 
parent haloes, the number of haloes with at least 10 satellites within r2oo, 
the number of satellite systems inclined by more than 80° to the galaxy and 
the number of these systems which have similar stellar masses M* (within 
a factor of 2) to the Milky Way and M3 1 (~ Mmw)- results are robust 
to more restrictive D/T cuts. Approximately 20% of satellite systems (with 
at least 10 members) are aligned perpendicularly to their pai'ent galaxy. 
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Figure 7. The cumulative probability distribution for the location of satel- 
lites relative to the host minor axis. Here, we consider all satellites within 
r2oo and are not restricted to those systems with large numbers of satellites. 
5 is defined from the minor axis (i.e. 5 = is polar alignment, S = 90° 
is planar alignment). The blue distribution is relative to the short axis of 
the galaxy (or disc) and the red distribution is relative to the short axis of 
the dark matter halo (defined for r ^ r2oo)- There is a preference for pla- 
nar alignment relative to the dark matter halo but there is a much weaker 
correlation with the galaxy. 

the previous calculation, as we consider all of the satellite gal axies 
and stack them together (c.f . lBrainerdll2005l : lYang et al.ll2006l) . We 
show the cumulative probability distribution of 5, which is defined 
as the angle between the short axis of the host and the position vec- 
tor of a satellite galaxy (i.e. (5 = is polar alignment and 5 = 90° 
is planar alignment). The blue and red lines are relative to the mi- 
nor axis of the galaxy and the minor axis of the dark matter halo 
respectively. An isotropic distribution is shown by the dotted line. 
The satellites exhibit a roughly planar alignment relative to the dark 
matter distribution, in agreement with the right hand panel of Fig.|6] 
where we consider only systems with a large number of satellites. 
The satellites have a relatively weaker alignment relative to the in- 
ner disc, although there is a slight bias towards planar ali gnment. 
It is in teres ting that we find qu alitatively similar results to iBrainerj 
( l2005h and I Yang et all ilOOdt) who, owing to the small number of 
satellite galaxies per parent halo, use stacked samples of satellites 
to generate a probability distribution of their orientations with re- 



spect to their hosts. However, a direct comparison is difficult due to 
the different halo selection criteria used by these authors. 

4.2 Rotational support 

It has been suggested that the sat ellites of the Milky Way lie in 
a rota tionally supported disc (e.g. iKroupa et al.ll2005l . iMetz et al] 
|200^ . We calculate the rotational velocity in the plane perpendic- 
ular to the short axis of the satellite distribution (Csat), {^4,') and 
find the velocity dispersion in this plane (cr^/). Typical values of 
these quantities are: (v^i) — 40kms^^ and a^i = 150kms^^.The 
left-hand panel of Fig. [8] shows a histogram of the ratio of this net 
rotational velocity to the velocity dispersion. Note that here we are 
restricted to satellite systems with ten or more members and focus 
on the ten brightest satellites within r2oo- The majority of satellite 
systems are not rotationally supported and their kinetic energy is 
dominated by internal motions. However, there are a small fraction 
(~ 9%) of systems that show substantial rotational support with 
\{V4,)\/a^, > 0.8. 

Further investigation into these rotationally supported sys- 
tems of satellites suggest that a significant fraction of the satel- 
lites are accreted in groups (i.e. from similar directions at the same 
time). We give an example in the middle panel of Fig. [8] This 
system of satellites defines a highly flattened plane at 2 = 
(Saat ~ 0.25) and shows evidence of substantial rotational sup- 
port with ~ 0.9. The radial trajectories of some of 
the satellites belonging to this system are shown in Fig. [8] We 
can see that 4 of the satellites infall as a group at 2 ~ 0.6 (or 
tiookback ~ 6Gyr). The coherency of the angu lar momenta of thi s 
group of satellites is retained until 2 = (c.f. iLi & Helmill2008h . 
Hence, signatures of rotational support are closely linked to sys- 
tems of satellites where a substantial number of the sate llite popu- 
lation today were accreted as a group. iMetz et al.l ( l2009h argue that 
the majority of the Milky Way satellites did not enter the halo in 
a group based on comparisons with dwarf associations in the Lo- 
cal Group (found by TuUv et al. 2006). Our result suggests a link 
between group infall and rotational support but does not directly 
addr ess the issue of gr oup infall into the Milky Way itself. 

iMetz et al.] j2008l) remark that at least 3 of the classical Milky 
Way satellites have orbital poles which lie within 30° of the short 
axis of the so-called disc of satellites. In the right-hand panel of 
Fig.jS] we show the distribution of the number of satellites with or- 
bital poles within 15°, 30° and 45° (red, black and blue lines) of 
the short axis of their spatial distribution. In ~ 20% of the systems, 
3 satellites have orbital poles within 30° of the normal to their spa- 
tially defined distribution. However, no systems have more than 7 
(out of 10) satellit es which all have orbita l poles within 30°. Thus, 
in agreement with lLibeskind et ah liooi), we find that an arrange- 
ment whereby the majority of the classical Milky Way satellites 
have orbital poles well aligned with the short axis of their flattened 
distribution will be difficult to reconcile with the results presented 
here. For comparison, we show with the dashed black line the case 
for an isotropic distribution of satellites in both positional and ve- 
locity space. The fraction of systems with more than two satellites 
with orbital poles within 30° of the short axis of their spatial dis- 
tribution is reduced in the isotropic case. This reinforces the point 
that the satellites are not randomly distributed in phase space. 

4.3 Angular Momentum Orientation 

The specific angular momenta of the satellites are computed us- 
ing eqn lO. In this case, we consider our whole sample of satellite 
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Figure 8. Left-hand panel: The probabihty distribution,P(| (V^/ ) |/(t^/ ), for the ratio of rotational versus dispersion support for each system of satellites 
where Ngi^tir ^ ?'20o) > 10. We only consider the ten brightest satellite galaxies within r200- The median rotation velocity {{V^')) and dispersion (a^i) 
are calculated in the plane perpendicular to the short axis of the satellite distribution. Middle panel: The radial trajectory of a system of satellites exhibiting 
rotational support at z = 0. A group of satellites is accreted at z ~ 0.6 (red-dashed lines). This group dominates the signature of rotational support shown at 
z = 0. Dotted lines show the trajectories of some of the other satellites and the solid blue line indicates r2oo(^). Right-hand panel: The fraction of systems 
{N / Ntata.l) where the orbital poles of Nsat satellites lie within 15° (red line), 30° (black line) and 45° (blue line) of the short axis defined by the spatial 
distribution of satellites. The horizontal black line indicates the fraction of systems which have 3 satellites with orbital poles within 6 < 30°. The dashed 
black line shows the case for 6 < 30° for an isotropic distribution of satellites (in both positional and velocity space). 



galaxies and are not restricted to parent haloes with a large number 
of satellite galaxies. 

The top panels of Fig. |9] show the distribution of the angles 
between the angular momentum vector of each individual satellite 
and the angular momentum of the dark matter halo (jdm {r ^ r2oo ) , 
left panel) and the angular momentum of the galaxy (jgai(T' ^ 
0.1r'2oo)i right panel). There is an obvious bias in both cases to- 
wards alignment between the angular momentum vectors. With 
respect to the angular momentum of the disc (dark matter) ap- 
proximately 61% (68%) of the satellites are on 'prograde' orbits 
(cos^ > 0) and 39% (32%) are on 'retro grade' orbits ( cos0 < 0). 
This is in agreement with recent work bv lLovell et al.l ( 1201 ih who 
analysed the orbital angular momentum of dark matter subhaloes in 
the Aquarius simulations. All six of the parent haloes in this study 
contain populations of co-rotating subhalo orbits. In addition three 
of their parent haloes contain subhaloes on retrograde orbits. The 
authors attribute these configurations to the filamentary accretion 
of subhaloes. 

Whilst we see a bias towards co-rotating orbits relative to the 
inner galaxy, the spatial orientation of the satellites is only weakly 
related to the orientation of the disc. For most satellite systems, 
there is little correlation between their shape and net angular mo- 
mentum (only a small fraction are rotationally supported). In addi- 
tion, phase space mixing is more rapid in spatial coordinates than 
velocity space — the angular momentum orientation of the satel- 
lites will be preserved for longer than any spatial anisotropy. This 
is especially true for satellites accreted at earlier times. 

In the bottom panels of Fig.[9l we show the radial dependence 
of the alignment between the satellite angular momenta and the 
dark matter halo angular momenta (left panel) and the galaxy an- 
gular momenta (right panel). The angular momenta of the satellite 
galaxies is more closely aligned with the dark matter than the inner 
galaxy. The continual accretion of material can alter the orientation 
of the net spin of the halo at larger radii. The accretion of satellites 
is closely linked to the spinning up of the dark matter halo (e.g. 



IVitvitska et aLll2002l) . hence it is unsurprising that the satellite an- 
gular momenta are well aligned with the dark matter spin. We found 
in Section 3.2 that there can be substantial misalignments between 
the angular momentum of the galaxy and the net spin of the outer 
dark matter halo. Hence, the angular momenta of the satellite galax- 
ies are less strongly related to the galaxy (this is especially true for 
those satellites accreted at later times). 

We find that there is a more pronounced bias towards satel- 
lites on co-rotating orbits (relative to both the disc and dark matter 
halo) in the inner regions of the halo. We relate this radial depen- 
dence to the accretio n history of the satellite galaxies in the follow- 
ing section. Note that lHwang & ParkI ( 1201 Ol) use a sample of SDSS 
galaxies to study the rotation of satellite galaxies relative to the spin 
of their host galaxy. The authors also find a stronger bias towards 
co-rotating orbits in the inner regions of their host galaxies. 



4.4 Anisotropic Accretion 

The redshift z — population of satellites is preferentially aligned 
in the plane perpendicular to the short axis of the dark matter distri- 
bution and exhibits a bias towards co-rotating orbits (with respect 
to both the net spin of the dark matter and the inner galaxy). Both 
these observations hint that the satellit e galaxies are accreted from 
preferential directions (e.g. filaments, iLibeskind et al.ll2005h . Ac- 
cording to tidal torque theory, the angular momentum of a halo is 
acquired through the tidal interactions between neighbouring struc- 
tures (out to 20Mpc) and infalling subhaloes. Presumably some of 
the surviving population of satellites today were accreted in a sim- 
ilar fashion and from similar directions as those subhaloes which 
originally spun up the halo. Thus, the bias towards co-rotating or- 
bits can be explained by tidal torque theory whereby the satellite 
population today bears the imprint of those subhaloes which were 
accreted from similar directions as the substructures involved in the 
early stages of galaxy formation. 

We trace back our z — population of satellite galaxies to 
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Figure 9. Top panels: Alignment between the angular momentum vector of each individual satellite galaxy (jaat.i) and its parent dark matter halo 
(idm(^ ^ ^20o)i Isft panel) and galaxy Ogai('" ^ 0.1r2oo) right panel). Downward pointing arrows give the median values and the eiTor bars denote 
Poisson uncertainties. Bottom panels: The median alignment values between the angular momentum vectors of the satellite galaxies and their parent dai'k 
matter halo (left panel) and and galaxy (right panel) as a function of radius. Error bars and shaded regions encompass the values covered by 68% of the 
distribution. There is a bias towards co-rotation (relative to both the spin axis of the dark matter halo and the spin axis of the disc) which is more pronounced 
in the inner regions of the halo. 



Our results are summarised in Fig.[TO] Note that we align 
our coordinate system with respect to the galaxy &l z = Q and jgai 
is defined at 2 = 0. The top left and middle panels show the per- 
centage of prograde (red) and retrograde (blue) satellite orbits as 
a function of the radial distance from the parent galaxy at z = 
and as a function of the time at which the satellite was accreted. 
We define the 'accretion' time when a satellite first crosses over 
'"200(2)- Similarly, the bottom left and middle panels show the me- 
dian misalignment between the galaxy angular momentum (defined 
at 2 = 0) and the satellites' angular momenta as a function of radial 
distance and accretion time. We see a stronger bias towards pro- 
grade orbits for those satellites accreted at larger lookback times. 
In addition, those satellites located closer in to the parent halo have 
a stronger tendency to be prograde that those located further out. 
The top right hand panel shows the median satellite Vmax, the peak 
of the satellite's circular velocity profile, at the time of accretion 
as a function of the radial distance from the parent galaxy today. 
Satellites located closer into the parent halo are on average more 



massive at the time of accretion. The bottom right hand panel shows 
the median alignment angle between the satellite angular momenta 
and the galaxy as a function of the satellite's Vmax at the time of 
accretion. The most massive satellites at accretion tend to have a 
stronger alignment with the net angular momentum of the galaxy. 

These results support the suggestion that the bias towards co- 
rotating orbits is a consequence of the hierarchical assembly of 
galaxies. We find that more massive satellites accreted at earlier 
times are more strongly biased towards co-rotation. These satel- 
lites resemble the early substructures which originally spun up the 
halo. The stronger bias towards co-rotating satellites in the inner 
regions of the halo can be explained by the prevalence of the most 
massive satellites at accretion in these regions at 2 = 0. By dynam- 
ical friction effects, we expect the more massive satellites to sink 
into the centre of the parent halo on shorter timescales than those 
with lower masses. 



^ We do not trace our sample of satellites beyond redshift 2 = 2 as only 
a very small fraction of satellites accreted before this epoch have survived 
until 2 = 0. 



In Fig. [TT] we show the 2D projection of the infall direction 
of the satellite population at 2 = 0. In the left-hand panel, 6 and <j> 
are the polar and equatorial axes defined with respect to the plane 
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Figure 10. Top left panel: The percentage of satellites on prograde (red) and retrograde (blue) orbits as a function of radial distance from the parent halo 
at 2 = 0. The dashed horizontal lines denote a ratio of 60:40. Error bars give Poisson uncertainties. Bottom left panel: The median alignment between the 
satellite angular momenta and the galaxy as a function of radial distance aX z = 0. The error bars indicate the dispersion in the distributions. Top middle 
panel: The percentage of satellites on prograde (red) and retrograde (blue) orbits as a function of accretion time. Bottom middle panel: The median alignment 
between the satellite angular momenta and the galaxy as a function of when they were accreted. Top right panel: The median satellite Vmax at the time of 
accretion as a function of radial distance from the host at z = 0. Bottom right panel: The median alignment between the satellite angular momenta and the 
galaxy as a function of satellite Vmax at the time of accretion. A stronger bias towards prograde orbits is seen by satellites accreted at early times and by those 
satellites which were relatively massive at the time of accretion. 
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Figure 11. The projection in cos(7r/2 — 6) and (j> of the direction of satellite accretion. 6 and ij> are the polar and equatorial axes as defined in the disc plane 
(left panel) at 2 = or defined relative to the dark matter halo shape (right panel) at z = 0. The right and bottom inset panels give the fractional number of 
satellites accreted as a function of 6 and as a function of (f>, respectively. Note that the dark matter halo shape is defined for r r2oo . The positions on the 
sky are calculated when a satellite passes inside r2oo (-2) for the first time. The scale bar gives the number of satellites in each pixel. There is a preference for 
satellite accretion in a plane perpendicular to the short axis of the dark matter halo (as defined at z = for r r200 ) but there is a very weak correlation with 
the galaxy. 
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Figure 12. The alignment of the axes of the satelhte velocity anisotropy with the galaxy (left panel) and the dark matter halo within r2oo (right panel). The 
shape axes of the inner galaxy are defined within the region r ^ 0.1r200 whilst the shape axes of the dark matter distribution are defined for r ^ r2oo- The 
velocity dispersion of the satellites tends to be maximum in the plane peipendicular to the short axis of both the dark matter halo (right panel) and the galaxy 
(left panel). 



of the discQ In the right-hand panel, 6 and cj) are defined relative 
to the dark matter shape (the shape is computed for r 7-200 )• We 
define a spherical surface at r — r2oo{z) and calculate 6 and (j> 
when a satellite crosses this surface. We stack all the haloes in our 
sample where 9 — defines either the plane of the disc (left panel), 
or the plane perpendicular to the short axis of the dark matter halo 
(right panel). There is a preference for accretion in a plane perpen- 
dicular to the short axis of the dark matter distribution. The bottom 
inset panel shows the fractional number of satellites accreted as a 
function of <^ - this shows there is a preference for accretion along 
the major axis of the dar k matter halo (i.e. \(f>\ = 0, 180° or ±x). 
This is in agreement with lLibeskind et al] ( l2005h and lZentner et al.l 
( l2005h who find a preferentia l direction of satellite accre tion along 
the major axis of the halo and lBailin & SteinmetjJlOOSbl) who find 
that there is a strong tendency for the minor axes of haloes to lie 
perpendicular to large-scale filaments. There is a much weaker bias 
relative to the disc plane. In Section 143] we found that the spatial 
orientation of the satellite galaxies at 2; = is much more strongly 
correlated with the dark matter shape rather than the galaxy. Note 
that both the disc plane and short axis of the dark matter halo are 
defined at z — 0. We note that this does not take into account the 
evolution of the orientation of the galaxy or dark matter halo. How- 
ever, our analysis does show that when satellites are accreted they 
have a preferential alignment relative to the dark matter halo as 
defined today. 



4.5 Velocity Anisotropy 

By restricting attention to those systems which have ten or more 
satellite galaxies, we now examine the properties of the velocity 
dispersion tensor The velocity anisotropy of actual satellite galaxy 
populations has never been measured, but it is an important param- 
eter for studies which use satellites as tracers of the dark matter 
potential. Fig[T2] shows the alignment of the principal axes of the 
velocity dispersion tensor with the shape axes of the galaxy and the 
dark matter halo. 



* The angle </) is defined relative to the major axis of the galaxy but note 
that as b/a ~ 1 for a disc this is a rather ambiguous definition 
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Figure 13. Top-left panel: The root-mean-square radial {Vr, red) and tan- 
gential (Vt, blue) velocity components as a function of radius. The solid and 
dashed lines are for the satellite galaxies and dark matter particles respec- 
tively. Top-right panel: The root-mean-square equatorial {V^ , magenta) and 
polar (Vg, green) components as a function of radius. Bottom-left panel: 
The velocity anisotropy parameter (fi) as a function of radius. /3 is consis- 
tent with isotropy over a wide radial range Bottom-right panel: The polar 
and equatorial velocity anisotropy components (Pg, fi^) as a function of 
radius. 



The short axis of the velocity anisotropy tensor tends to align 
with the short axis of the galaxy (and the short axis of the dark 
matter distribution). This suggests the velocity dispersion of the 
satellites tends to be maximum in the plane perpendicular to the 
short axis, i.e. in the plane of the disc. This cylindrical alignment 
of the velocity anisotropy tensor agrees with our earlier findings 
that there is a bias toward co-rotating satellite orbits. Note that here 
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Figure 14. The polar (/Sq, right panel) and equatorial (/3^,left panel) velocity anisotropy parameters against the overall anisotropy parameter, /3. Each point 
represents a system of ten or more satellites (within 2r2()o) belonging to an individual galaxy. Red, blue and green points represent orbits dominated by radial, 
equatorial and polar motions respectively. There is significant halo-to-halo scatter in the overall velocity anisotropy. 



we have no direct comparison to the Milky Way satellites as their 
tangential velocity components are poorly constrained. 
The anisotropy parameter, /3 is defined as 



/3= 1- 



= 1 - 



{vj} + jv') 



(9) 



This is a measure of how radially or tangentially biased the satel- 
lite orbits are. To distinguish between polar and equatorial biased 
orbits, we define two more anisotropy parameters 



{vi) 



P4. 



(10) 



These quantities are defined in a cylindrical polar coordinate sys- 
tem {R, (j), z) aligned such that the z-axis is normal to the disc. 

In the top-left panel of Fig.[T3] we show the radial dependence 
of the radial and tangential velocity components. The solid lines 
are for the satellite galaxies, whilst the dashed lines are for the dark 
matter. Radial velocities dominate over tangential velocities at all 
radii for the dark matter particles. The satellite galaxies have more 
tangentially biased orbits in the inner regions (r < 0.5r2oo) and 
only become radially biased at larger radii (r > r2oo)- Over a wide 
radial range, the satellites are consistent with an isotropic veloc- 
ity distribution. This can also be seen in the radial dependence of 
/3 shown in the bottom-left panel. The two right-hand panels de- 
compose the tangential velocity component into its polar and equa- 
torial parts. The tangential bias of the satellite orbits in the inner 
regions are dominated by their equatorial motion. This is seen in 
the bottom-right panel where is (comparatively) large and neg- 
ative at small radii. The equatorial component dominates over the 
polar component for the dark matter particles over the whole ra- 
dial range. Dark matter haloes are dispersion supported (see Fig. 
[2l(, and the dominance of equatorial velocity components over the 
polar components leads to the (slight) flattening of the dark matter 
halo in the z-directior[3. 



^ Note the flattening of the dark matter halo is not always in the z-direction 
(defined relative to the inner disc) as there can be misalignments between 
the short axis of the dark matter distiibution and the z-axis of the galaxy, 
especially at larger radii (see Section[3). However, there is reasonably good 
alignment for the majority of haloes. 



In Fig. [TS] we have stacked all the satellites in our sample 
together by normalising radial distances by r2oo and velocities by 
V20()- For parent haloes with ten or more satellites within 2r2oo, we 
plot the velocity anisotropy parameters, /3, [ig and 13^ in Fig. [14] 
It is apparent that there is a wide spread in the /? parameters for 
individual haloes ranging from —2 < /? < 0.8. By stacking all the 
satellites of different haloes together, we find /3 values ranging from 
—0.4 < /3 < 0.2 over the same radial range. We can see from Fig. 
1 141 that there is a slight bias towards satellite systems dominated 
by motion in the equatorial plane (39%) but many satellite systems 
are dominated by their radial motions (35%) or are biased towards 
polar orbits (26%). 

The velocity anisotropy parameter, f5, is an important factor 
required to estimate the masses of local group galaxies such as the 
Milky Way and M31 (see next Section |5]l. This parameter is dif- 
ficult to observe as in most cases only line of sight velocities are 
available and we lack full 3D velocity information. To overcome 
this, many a uthors adopt anisotropy parameters d erived from simu- 
lations (e.g. lXue et a"Ll(2008l ; IWatldns et al.ll2010l) . Our finding that 
there is significant halo-to-halo scatter means that significant cau- 
tion is warranted when applying a velocity anisotropy applicable to 
a simulated satellite system of an individual halo or from a stacked 
sample of satellites to our own Milky Way Galaxy. 



5 AN APPLICATION: MASS ESTIMATORS 

In practice, simple estimators are often used to compute the mass of 
a dark halo from the positions and velocities of the satellite galax- 
ies (e.g., Watkins et al. 2010). They depend on simplifying assump- 
tions, such as underlying spherical symmetry of the dark halo, or 
constant velocity anisotropy. As we have seen, there are numerous 
effects present in the simulations - triaxiality of the halo, contin- 
uing infall to the present day, variation of anisotropy with radius 
- that are not accounted for in the mass estimators. Hence, it is 
interesting to see how the estimators fare against simulation data. 

The projected m ass estimator (hereafter PME) introduced by 
iBahcall & Tremaind ( IT98Ii) takes the form; 



(11) 
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Figure 15. Histograms showing the ratio of the estimated mass to the true 
mass for the two different mass estimators. The projected mass estimator 
(PME) and Tracer mass estimator (TME) are shown by the blue and red 
lines and respectively. We give the fraction of reasonable estimates (FRE) 
and the inter-quartile range (IQR) as defined in the text. 



where R and vios are the projected positions and line-of-sight ve- 
locities of the A'^ satellite galaxies. The constant is Cp — 16/tt 
(isotropic) or Cp = 32 /n (radial orbits). 

The tracer mass estimator (hereafter TME) is given in Watkins 
et al. (2010, see also Evans et al. 2003). It assumes a spherically 
symmetric power-law for the halo potential $ oc r~". We use the 
form: 



M 



Ct 
GN 



E 



(12) 



where a ~ 0.5 and the constant Ct is given in eqn (26) of Watkins 
et al. (2010) and depends on the velocity anisotropy /3. In particu- 
lar, $ oc r~°'^ is a good approximation to a Na varro-Frenk- White 
(NFW lNavarro et alJ[l99a : lNavarro et alJll997l) profile often used 
to describe simulated dark matter haloes. We assume isotropic or- 
bits (/? = 0), but investigate the validity of this assumption by com- 
paring the estimated masses when the true velocity anisotropy of 
the tracer satellites is used. 

Eqns dl lb and il2\ provide an estimate for the total mass 
within the radius of the furthest tracer (rout). We select an arbitrary 
viewing angle for our simulated haloes to generate projected posi- 
tions and line of sight velocities. We then compute the 'true' mass 
within r-Qut for each halo and compare to masses found via the two 
mass estimators. Note we only select satellites within r < 1.5r2oo- 
We use all satellites, but check that our results are not significantly 
affected when only luminous satellites are included 

In Fig. [15] we present histograms of the ratio between the es- 
timated and true mass. We define the 'Fraction of Reasonable Es- 
timates' or FRE as the fraction of estimates in the range 0.5 < 
Mest/Mtruo < 2. We also give the Inter Quartile Range (IQR) 
of the mass estimates in Figure [15] Both the TME and the PME 
perform well, with ~ 81% and 78% of the estimates satisfying 



our 'reasonable' criteria, respectively|f[ The TME performs slightly 
better than the PME (note the more symmetrical distribution in 
Fig. |15l l, but both provide good results especially considering the 
rather idealized assumptions under which they are derived. The 
IQR shows that the uncertainty in the mass estimates given by the 
simulations is of si milar magnitude (~ 30%) to the statistical un- 
certainty found by I Watkins et al.l ( 1201 Ol) in their estimates of the 
Milky Way mass. 

Fig.[T6]illustrates how the mass given by the TME varies when 
we use the true velocity anisotropy parameter instead of assuming 
an isotropic distribution. Assuming isotropy for tangential orbits 
leads to an overestimate in the mass, whilst the reverse is true for 
radial orbits. The median absolute difference between the isotropic 
estimator and the true anisotropy estimator is 2%Mtruo- Unless 
the velocity anisotropy is strongly radial or tangential, the assump- 
tion of isotropy yields reasonable mass estimates. This is impor- 
tant for the application of such estimators to observational data, 
as we rarely have observationally derived values of the velocity 
anisotropy parameter. 

In our samples of satellite galaxies, approximately 3% are un- 
bound. Whilst unbound satellites are not common, their inclusion 
in mass estimators can cause fairly large deviations from the true 
mass value. We compute the estimated mass with and without the 
unbound satellites for the 46 haloes which contain at least one un- 
bound satellite (within 1.5r2oo)- The left hand panel of Fig. [TT] 
shows the distribution of the relative difference between the mass 
estimates. The right hand panel plots the two mass estimates against 
one another. Inclusion of unbound satellites causes a systematic 
overestimate of the mass. The TME and PME have median differ- 
ences of 45%A/true and 55%Mtruc when unbound satellites are 
included. 



6 CONCLUSIONS 

We investigated the orbital properties of the sa tellites of late-typ e 
galaxies using the GIMIC suite of simulations l lCrain et alj[2009h . 
These state-of-the-art simulations incorporate baryonic physics 
into a ACDM cosmological framework and produce realistic disc 
galaxies at 2: = 0. We analyse the phase space distributions of 
the satellite galaxies relative to the luminous baryonic material (i.e. 
the central galaxy disc) as well as the unseen dark matter compo- 
nent with a large sample of galaxies. In this way we can provide 
a more direct comparison with observations. Our sample of par- 
ent haloes were chosen to be relaxed systems in the mass range 
5 X 10^^ < A/200 /A'/© < 5 X 10^^, which broadly overlaps with 
the mass of our own Milky Way galaxy and of M3 1 . 

The parent dark matter haloes in our sample are generally tri- 
axial but roughly spherical and have axial ratios which are roughly 
constant with radius ((sdm) ~ 0.8). Comparison to the dark matter 
only counterparts of our sample shows that the inclusion of bary- 
onic physics affects the shapes of the dark matter haloes signifi- 
cantly, even out to r-2oo- The central galaxy (or 'disc') is often mis- 
aligned in both shape and angular momentum with the dark matter 
halo. The inner regions of the dark matter halo are well aligned 
with the central galaxy but there can be substantial misalignments 
at larger radii, in the region most relevant for the satellite galaxies. 
We find that for radii r ~ r2oo, the short axis of approximately 



^ By restricting ourselves to systems with 10 or more satellites the fraction 
of reasonable estimates increases to ~ 90%. 
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Figure 16. Left hand panel: Histogram showing the absolute difference in the estimated mass when the actual velocity anisotropy is used instead of assuming 
/3 = 0. This is normalised by the true mass. Right hand panel: The variation in the mass estimator when the true velocity anisotropy is used as a function of 
velocity anisotropy. 




Figure 17. Left hand panel: Number distribution showing the relative difference in the mass estimates when unbound satellites are included or excluded in the 
analysis. Median values are shown by downward pointing arrows. Right hand panel: The ratio between the estimated mass and true mass in the case where 
unbound satellites are included (y-axis) and when they are not (x-axis). The colour scheme is the same as Fig. 1151 Arrows represent the median values. 



30% of our parent dark matter haloes are significantly misaligned 
from the short axis (or z-direction) of the inner galaxy (6 > 45°). 
In velocity space, the net spin of the dark matter halo can be almost 
perpendicular to the angular momentum vector of the inner galaxy 
in ~ 40% of our sample and 2% are even spinning in the opposite 
sense to the inner galaxy. 

There is an obvious spatial bias between the dark matter of 
the parent halo and the satellite galaxies. The satellite system has 
a more flattened shape than the dark matter halo and is not as cen- 
trally concentrated. By considering all systems with 10 or more 
satellites within r2oo , we find that the satellites preferentially align 
in a plane perpendicular to the short axis of the dark matter halo. 
However, owing to the misalignments between the inner galaxy and 
(outer) dark matter halo, this pref erential alignment is much we aker 
relative to the inner galaxy (cf. lAgustsson & Brainerj|200^ . In 
fact, by only considering the 10 highest stellar mass satellites in 
each system, we find the distribution of satellites is almost uni- 
form relative to the central disc and the 'unusual' orientation of the 
classical Milky Way dwarfs is not uncommon: 20% of satellite sys- 
tems are perpendicular (within 10°) to the disc. In a similar fashion 



to iBrainerj ( |2005|) . we find the probability distribution of the ori- 
entation of the satellite galaxies relative to their hosts by stacking 
all of the satellites in our sample. In qualitative agreement with 
the observational results from SDSS and the 2dF Galaxy Redshift 
Survey, we find that the satellites have a weak bias towards planar 
alignment relative to the disc. However, there is a much stronger 
alignment relative to the dark matter halo shape. 

It has been suggested that the Milky Way sate llites may oc- 
cupy a rotationally supported disc jMetz et al]|2008h . We find that 
satellite systems which are planar and rotationally supported are 
relatively uncommon (~ 9%) in the simulations. This often occurs 
when a large fraction of the satellites is accreted in a group that re- 
tains its coherence in velocity space until z — 0. We find that it is 
not unusual to have 3 out of 10 satellites with orbital poles within 
30° of the normal to their spatial configuration. This is consistent 
with the available proper motion data on the classical dwarfs. How- 
ever, we find that if a substantial number of the classical dwarfs of 
the Milky Way (e.g. 7 out of 10) were found to have orbital poles 
aligned with the normal to the disc of satellites, then this would be 
inconsistent with the results of our simulations. 
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We tested this claim in the simulations and found that satellite 
systems which are planar and rotationally supported are relatively 
uncommon (~ 9%). This often occurs when a large fraction of the 
satellites is accreted in a group that retains its coherence in velocity 
space until z = 0. We find that it is not unusual to have 3 out of 10 
satellites with orbital poles within 30° of the normal to their spatial 
configuration. This is consistent with the available proper motion 
data on the classical dwarfs. However, we find that if a substantial 
number of the classical dwarfs of the Milky Way (e.g. 7 out of 10) 
were found to have orbital poles aligned with the normal to the disc 
of satellites, then this would be inconsistent with the results of our 
simulations. 

There is a bias towards co-rotating satellite orbits relative to 
both the angular momentum of the disc and the net spin of the dark 
matter halo. This is more pronounced in the inner regions of the 
halo. This co nfirms earlier resu lts relating to dark matter only sim- 
ulations (e.g. lLovell et al l|), but the bias with respect to the in- 
ner disc is weaker owing to the angular momentum misalignments 
between the inner galaxy and dark matter halo. A preference for co- 
rotating orbits is a natural consequence of the hierarchical assembly 
of galaxies whereby satellites accreted at earlier times are related 
to those substructures that helped spin up the galaxy. We confirmed 
this by finding a stronger bias towards prograde orbits by the more 
massive satellites that were accreted at earlier times. By tracing 
back the infall orientation of our sample of satellite galaxies, we 
find their anisotropic distribution is due to their preferential accre- 
tion in directions perpendicular to the short axis of the dark matter 
distribution. There is a weaker correlation with the orientation of 
the inner galaxy. 

The velocity anisotropy tensor for the satellite galaxy systems 
is cylindrically aligned relative to the central disc. Tangential mo- 
tions dominate at smaller radii, often due to the prevalence of equa- 
torial (as opposed to polar) orbits. The velocity anisotropy, /3, is 
an important parameter which is largely inaccessible with present 
observations. Here we show that /3 is consistent with zero over a 
large radial range when all satellites are stacked together. How- 
ever, inspection on a halo-by-halo basis shows that there is a signif- 
icant degree of scatter between haloes. This scatter puts into ques- 
tion the validity of using a single simulation as a template velocity 
anisotropy to be applied to real galaxies. 

Finally, we considered an application of the orbital properties 
of the satellite galaxies. We tested two popular mass estimators in 
the literature which make use of the projected positions and line of 
sight velocities of tracers, such as satellite galaxies, to estimate the 
mass of the parent halo. The projected mass estimator (PME) and 
the tracer mass estimator (TME) both perform well and estimate 
'reasonable' (within a factor of 2 of the true mass) halo masses. 
The TME performs slightly better as it assumes the satellites are 
tracers of the halo potential rather than having their density profile 
generated by the dark matter potential (as assumed by the PME). 
We found that an unknown velocity anisotropy parameter can lead 
to incorrect mass estimates but these are only substantial when 
/3 is significantly non-isotropic. In addition, the inclusion of un- 
bound satellites can cause large overestimates of the true halo mass. 

Requests for simulation data should be directed to 
ajd75@ast.cam.ac.uk. 
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APPENDIX A: CONVERGENCE TESTS 

Our sample of simulated galaxies is drawn from the five 
intermediate-resolution GIMIC simulations, which have been run to 
2 = 0. In addition, there is one higher resolution GIMIC simulation 
(the -2a region) available to 2 = 0. This simulation has eight times 
better resolution and allows us to assess the numerical convergence 
of our results. We have checked that the main results in this paper 
are unchanged in this higher resolution GIMIC simulation. Here, we 
give examples for two of our main results. 

In Rig. lAll we we show the distributions of misalignments be- 
tween the galaxy and dark matter halo short axes both for r ^ 
0.1''2oo (left panel) and for r ^ r2oo (right panel) for the high res- 
olution gimic runs. This is shown in the bottom panels of Rig. [5] 
in the main text. There is very good agreement between the high 
and intermediate resolution simulations. Our finding in the main 
text that there can be significant misalignments between the galaxy 
and the outer dark matter halo is therefore robust to an increase in 
resolution. 

In Rig. IA2I we show the distribution of aligimients between 
the short axes of the satellite systems and the short axes of the dark 
matter halo (for r ^ r20()) and the galaxy of their parent haloes 
(see Fig. |6]in the main text). We only consider the 10 brightest 
satellites within r2oo- The black and blue lines show the distribu- 
tions in the intermediate and high resolution simulations respec- 
tively. Note that only the — 2cr volume is available at z = in the 
high resolution simulations so we do not achieve the same statis- 
tics as the intermediate resolutions runs. Approximately 30 parent 
haloes have 10 or more satellites within r2oo in the high resolution 
sample (c.f. ~ 80 in the intermediate resolution sample). Ror both 
the intermediate and high resolution runs the satellite distribution 
preferentially aligns in a plane perpendicular to the short axis of 
the dark matter distribution and the satellites show no preferential 
alignment relative to the galaxy. Thus, our conclusions from Rig.|6] 
in the main text are robust to an increased resolution. 
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Figure Al. The distribution of misalignment angles between the short axis of the galaxy and the short axis of the dark matter halo for r ^ 0.1r2oo (left panel) 
and r ^ r200 (right panel) respectively. This is for the high resolution GIMIC simulations (c.f. Fig.|3]for the intermediate resolution version of this plot). 




Figure A2. The orientation of the short axes of the satellite distribution relative to the short axes of the parent galaxy and dark matter halo (defined within 
?'20o)- The 10 brightest satellites within r200 are used to compute the shapes of the satellite distribution. Downward pointing arrows denote the median of the 
distributions and the dotted lines indicate a uniform distribution. The error bars denote Poisson uncertainties. The black and blue lines show the distributions 
in the intermediate and high resolution simulations respectively. 
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